A simple model for magnetic reconnection heated corona 
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ABSTRACT 

We construct a simple model for a magnetic reconnection heated corona above a thin accretion 
disk in AGNs and Galactic black hole candidates(GBHCs). The basic assumptions are that (1) the 
magnetic reconnection heat is cooled down overwhelmingly by Compton scattering in the corona 
and that (2) thermal conduction is dominantly cooled by evaporation of the chromospheric plasma 
in the disk-corona interface before Compton cooling sets in. With these two basic equations 
as well as equi-partition of magnetic energy with gas energy in the disk, we can consistently 
determine the fraction of accretion energy dissipated in the corona without free parameters, and 
thus determine the temperature and all other quantities in both the corona and disk for given 
black hole mass and accretion rate. Then, we calculate the luminosity contributed from the disk 
and corona and the coronal flux weighted Compton y parameter. It is found that, at a low 
luminosity (L < O.lLEdd) the spectrum is hard with spectral index of a sp ~ 1 (/„ oc ^~ Q5p ), 
while at a high luminosity (L > O.lLEdd) the spectrum can be either soft or hard. We also find 
that the situation is almost same for supermassive and stellar-mass black holes. These features 
are consistent with observations of AGNs and GBHCs. 

Subject headings: accretion, accretion disks-magnetic fields-galaxies: nuclei-X-rays:galaxies 



1. Introduction 

There is a common belief that the hard X- 
ray radiation of GBHCs and active galactic nu- 
clei (AGNs) arises from hot gas around black 
holes, either in forms of hot accretion flow or 
accretion-disk corona (e.g. Liang & Nolan 1984; 
Mushotzsky, Done, & Pounds 1993; Narayan, Ma- 
hadevan, & Quataert 1998). The structure of 
hot accretion flow has been extensively investi- 
gated recently based on the vertically one-zone ap- 
proximation. Many calculations have so far been 
presented and various notions, such as ADAF, 
CDAF, and ADIOS, have been proposed (Ichi- 
maru 1977; Blandford & Begelman 1999; Igumen- 
shchev & Abramowicz 2000). Such studies can be 
directly compared with observations. The ADAF 



model, for example, can reasonably fit the ob- 
served spectra of low luminosity AGNs and X- 
ray novae in quiescence. Nevertheless, co-existing 
soft and hard X-ray components observed in ordi- 
nary AGNs and at least some GBHCs cannot be 
well described by simple one-zone hot flow mod- 
els and require composite disk-corona structure. 
Yet, we do not have a standard way of model- 
ing disk corona despite many attractive attempts 
(e.g. Haardt & Maraschi 1991; Meyer, Liu & 
Meyer-Hofmeister 2000; Rdzahska & Czerny 2000; 
Kawaguchi, Shimura, & Mineshige 2001). One of 
the reasons is that, in addition to complex radia- 
tion and energy interaction between the disk main 
body and corona, magnetic fields, which tend to 
produce time variations and spatial inhomogene- 
ity (Kawaguchi et al. 2000), seem to play essential 
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roles there and are thus difficult to model in a sim- 
ple way. 

There are good reasons to believe that accre- 
tion disk corona may share some common fea- 
tures with solar corona (e.g., Galeev, Rosner, & 
Vaiana 1979; Zhang et al. 2000). For example, 
violent fluctuations are observed in both coronae 
(Takahara f 979) and the fluctuation properties are 
remarkably similar in the power-spectral density 
plot (UeNo et al. f997). Interesting theoretical in- 
vestigations have also been presented in the frame 
of magnetized disk corona (e.g. Di Matteo 1998; 
Di Matteo, Celotti, & Fabian 1999; Miller & Stone 
2000; Merloni & Fabian 2002). Then, to what 
extent the accretion disk corona is similar to the 
solar corona in terms of basic physical processes? 
What can we learn from the study of the solar 
corona for modeling the accretion disk corona? 

There has been a long history of the solar 
corona research, and recently our knowledge on 
its structure and time variations has been enor- 
mously enriched thanks to the intensive X-ray ob- 
servations from the space by the Yohkoh satellite 
(Uchida, Kosugi, & Hudson 1996). To understand 
the observed emission measure and temperature 
relation (EM-T relation), Shibata & Yokoyama 
(1999) propose a magnetic reconnection model for 
solar flares. They point out that only two equa- 
tions suffice to explain the observed EM-T rela- 
tion, which are the energy balance between recon- 
nection heating of magnetic field and conductive 
cooling, and the energy equi-partition between gas 
and magnetic fields. They find the theoretical EM- 
T relation in good agreement with the observation 
of solar and stellar flares over more than 10 orders 
of magnitudes in EM (see also Yokoyama & Shi- 
bata 2001). Furthermore, the coronal density can 
be estimated from the balance between conductive 
heating and evaporation cooling of chromospheric 
plasma, which is also consistent with observations. 

With this potentially powerful tool of describ- 
ing the basic properties of the solar corona, we are 
tempted to understand the accretion disk corona 
along a similar line. This is the starting point of 
this Letter. Our aim is to propose a simple and 
self-consistent method, based on the interaction 
between the disk and corona, to describe the basic 
properties of the disk and corona for given black 
hole mass and accretion rate, with no free parame- 
ters constrained by observations. This is a big dif- 



ference of our attempt from previous magnetized 
disk corona models. 

2. Magnetic reconnection heated corona 

First, we consider what relations determine the 
thermal structure of magnetic reconnection heated 
corona. From the X-ray observations of GBHCs 
and AGNs we know that the emitting temper- 
ature of hot plasma is ~ 10 9 K, and the power- 
law spectrum requires a Compton y parameter of 
y = -^^T ~ 1, leading to the constraint on the 
density of hot plasma to be n < 5 x 10 9 cm~ 3 in a 
hot region of size I ~ 10i?s {Rs is Schwarzschild 
radius). A successful model should account for 
these observational data. 

The energy balance in the magnetic flux tube 
holds between heating by magnetic reconnec- 
tion and cooling by Compton scattering, ther- 
mal conduction, cycle-synchrotron radiation. For 
a typical AGN system with a black hole mass 
of 10 s Mq and an accretion rate of O.lMEdd 
(-^Edd = lO^Edd/c 2 ), the radiation energy den- 
sity in disk is U la d « 10 5 ergs cm" 3 . The magnetic 
filed is B ~ 10 3 G under the assumption of energy 
equi-partition with gas in the disk. We find that 
the Compton scattering is the most efficient cool- 
ing mechanism in the range of above parameters 
(which is also true in GBHCs shown in Sect. 4). 
We thus have 

B 2 4kT 

-t~Va ~ ^na T cU ra ,d£, (1) 

47r m e c 

where Va is the Alfven speed, J7 ra d = U™ d + U™ d 
the soft photon filed to be Compton scattered 
from both the intrinsic disk and reprocessed ra- 
diation, I the length of magnetic loop, and other 
constants have their standard meanings. There- 
fore, the temperature in the magnetic flux tube is 
T = 1.21 x l^n^ ,2 U^Blt^K, which is around 
the observed value of <~ 10 9 K for typical param- 
eters of rag, U5, B3, and £14 in units of 10 9 cm~ 3 , 
10 5 ergs cm~ 3 , 3 x 10 14 cm(= 10i?s), respectively. 

Then how to determine the density? At the 
bottom of the flux tube, structural calculations 
for both solar corona (Shmeleva & Syrovatskii 
1973) and disk corona (Liu, Meyer, & Meyer- 
Hofmeister 1995; Meyer et al. 2000, Liu et al. 
2002) show that there is a very thin transition 
layer between the corona and the chromosphere, 
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where temperature changes very steeply within 
a depth h (h <C £)■ The thermal conductive 
flux there is quite large. If the initial mass den- 
sity is too low, the conduction rate can exceed 
the Compton cooling rate within this thin layer 
(i.e. conduction time scale is less than Compton 
timcscale). This conduction flux then heats up 
some of the chromospheric plasma into the mag- 
netic tube in a similar way to the solar corona 
(Yokoyama & Shibata 2001). This process is called 
chromospheric evaporation. The density of evapo- 
rated plasma is estimated from the energy balance 

1 /2 

at the interface, « ^n cvap /cT(^) 

(where the evaporating speed is around sound 
speed due to the large pressure gradient, k ~ 
lO^ergcm-^-^- 7 / 2 , 7 = 5/3). This plasma 
eventually diffuses to the whole loop, and the 
averaged density in the magnetic tube n fulfills 
n cvap h = nl. We thus get, 



k Ti 



7 



M-) 

7-1 \m H J 



(2) 



Once the density n reaches a certain value, the 
Compton cooling sets in and finally overwhelms 
the evaporation cooling, an equilibrium is then 
established between the magnetic heating and 
Compton cooling through the whole loop, evapo- 
ration at the interface becomes quite slow or even 
stops. For a further check of the argument, we 
compare the enthalpy flux with conductive flux 
from numerical calculations on the vertical struc- 
ture of the frictionally heated disk corona (Liu et 
al. 2002), finding that the approximation is cor- 
rect in order of magnitudes. For T = 1.21 x 10 9 K, 
the density n given by Eq.(2) is 1.6 x 10 9 cm~ 3 , 
which is just what is required. This implies that 
the evaporation can build up a corona with a den- 
sity ~ 10 9 cm~ 3 . 

Therefore, we expect two phases caused by 
the magnetic reconnection heating: firstly, mass 
evaporation at the bottom of magnetic flux tube 
quickly builds the corona up to a certain den- 
sity; then, Compton scattering in the tube steadily 
radiates away the magnetic heating. Combining 
Eq.(l) with Eq.(2), we derive T and n in the 
corona, 



t = 1.02 x io 9 c/; 
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n = 1.12 x 10 y f7 5 2 B 3 2 £ 14 4 cm- 



(3) 
(4) 



Eqs.(3) and (4) reproduce the typical temperature 
of 10 9 K and a density of 10 9 cm~ 3 and establish the 
relations between the disk and the corona thor- 
ough disk radiation and magnetic field. 

The magnetic energy is assumed to be equi- 
partitioned with gas energy of the disk, [3 = 
fidisk&^disk/ = 1. Other possible values of (3 
have been studied by Miller & Stone (2000). In 
our model, larger (3 means less energy dissipated 
through the corona, and results in lower coronal 
T and n as shown in Eqs.(5)-(8). To calculate 
disk quantities the classical disk model (Shakura & 
Sunyaev 1973) is used, in which absorption is dom- 
inated by electron scattering and pressure by ei- 
ther gas pressure f j\ sk or radiation pressure -PJ isk , 
depending on the accretion rate. For Pdisk = -PJisk' 
we find that the soft photon field is dominated by 
the intrinsic disk radiation, and T and n in the 



corona are 



T = 1.36 x 10 9 a 3 32 ft s m~ 



0.1' 10 *-io 



£? Q K (5) 



and 



2.01 x 10 9 a 



0.3 



'ft" 



0.1' 10^10** Cm 



-3 



(6) 

with standard parameters, viscous coefficient ao.3, 
ft, mg, fho.i, r io, ho m units of 0.3, 1, 10 8 M Q , 
0.1M E dd,10i?s, 10-Rs, respectively. For P disk = 
P| isk , the soft photon field is dominated by the 
reprocessed coronal radiations [/ ra d w OAUb (~ 
40% of the coronal radiation is reprocessed, see 
Haardt & Maraschi 1991) and we have 

T = 4.86 x lO^fft^mf m^r^ l\ K , (7) 

n = 2.55 x lO^aofft^^m^Th^r^^cm- 3 . 

(8) 



3. Fraction of accretion energy released 
through corona 

So far, we have not considered a back reaction, 
that is, the energy transfered from the disk to the 
corona by magnetic reconnection essentially af- 
fects the disk structure. Defining / as the fraction 
of accretion energy released in the reconnected 
magnetic corona, 

f _F COT (B\\ ( ZGMM Y 1 

f= F- = {^ VA ){-^Wr) ' (9) 
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we find that, in the expression of disk quantities, 
such as pressure (and hence B) and f7 ra d, ^o.i 
should be replaced by (1 — /)mo.i. Then the coro- 
nal temperature and density in Eqs.(5)-(8), should 
be modified appropriately, which, in turn, affects 
the amount of reprocessed radiation. Finally, we 
obtain an equation concerning / from Eq.(9) for a 
radiation pressure-dominant disk 

„ _M _£ _JL o 225 3 

/ = 1.45(1 - /)" 2 a .f ft s m 8 "ro^rtf ff , 

(10) 

and for a gas pressure-dominant disk 

. n _99 _11 11 J_ _JLL 3 

/ = 4.70xl0 4 (l-/)TOa |Oft « ^80^10^^160^ 

(11) 

Therefore, the ratio / is not a free parameter in 
our model and can be calculated out by solving 
Eq.(lO) or Eq.(ll) for given black hole mass and 
accretion rate. 

With < / < 1, Eq.(lO) has either double 
solutions or no solution. Only at small distances, 
R < Riad, there exist double solutions, where 

i?rad = SRsal^mprh^e'* . (12) 

This is to say, the radiation pressure-dominant 
disk exists only in the innermost region at rel- 
atively high m. At low to, the critical radius 
[Eq.(12)] becomes less than 3Rs, the radiation 
pressure-dominant region would no longer exist. 
In contrast, Eq.(ll) always has unique solution, 
indicating that the gas pressure-dominant region 
can extend to very small radius even at high to as 
long as a large fraction of accretion energy dissi- 
pates in the corona (/ <~ 1). 

Fig.l shows the radial distribution of variables 
T, n, B, V A , and r for M = 0.1M Edd (left) and 
M = 0.5M E dd (right). At 0.1M Edd there is no ra- 
diation pressure-dominant disk even close to the 
last stable orbit, and almost all the energy is trans- 
fcrcd to the corona and radiated away by Comp- 
ton scattering. The temperature in the corona 
is - 2 x 10 9 K, the density - 3 x 10 9 cm" 3 , the 
magnetic field ~ 10 3 G, optical depth ~ 0.8, and 
Alfven speed ~ 0.2c. The magnetic pressure be- 
fore reconnection is ~ 30 gas pressure of the recon- 
nected corona, which indicates a ~ 30 times verti- 
cal decrease of gas pressure from the disk midplanc 
to the corona under the equi-partition principle. 
With increase of M, all these quantities increase 



slightly until the radiation pressure-dominant re- 
gion appears in the innermost region. Temper- 
ature and other quantities in this state are lower 
than that at gas pressure-dominant state as shown 
in right panel of Fig. 1 . 

4. Compton y parameter and observa- 
tional implications 

The shape of the Compton spectrum depends 
on the y parameter. In an optically thin corona 
with thermal-distribution non-relativistic elec- 
trons, we have y = ^ fc ^ 2 naxi and spectral in- 
dex a sp = yf + f — §• Fig-1 shows that the 
Compton y parameter slowly decreases outwards 
along distance. In order to compare the spec- 
trum with observations we calculate the coronal 
flux weighted y parameter from 3i?s to 50i?s for 
a series of mass accretion rates. Mean fraction / 
(weighted by total accretion energy) of accretion 
energy radiating from the corona is also calcu- 
lated so that we can know the relative strength 
of luminosity from disk and corona. Results plot- 
ted in Fig.2 show that y -0.7 at M = 0.1M Edd , 
and that the accretion energy is overwhelmingly 
dissipated in the corona (/ « 1) with a coronal 
luminosity L COI w 4.3 x 10 44 ergss _1 . At a rate 
of 0.5M Edd , radiation pressure-dominant region 
appears in the innermost part, where / « 0.23, 
y w 0.41, and L COI w 5.1 x 10 44 ergss _1 . At an 
even higher accretion rate, M — M Ec id, the ra- 
diation pressure-dominant part expands and the 
emission of the system is dominated by the disk. 
We find / fts 0.03 and the coronal luminosity is 
mainly contributed from outer region (40i?s to 
50i?s) , i cor ~ 1-5 x lO^ergss" 1 , which is smaller 
than that at M = 0.1M E dd, and the correspond- 
ing 0.23) is lower. The results imply that, 
for a high rate the emergent spectrum is soft: a 
large fraction of the accretion energy is dissipated 
in the disk, the system is bright in optical-UV and 
soft X-rays contributed by the strong disk, but 
is faint in hard X-rays contributed by the weak 
corona, and thus the spectrum is steep. While for 
a low mass accretion rate, the spectrum is hard: 
almost all the accretion energy is channeled to the 
magnetized corona and is radiated by Compton 
scattering, the y parameter is 0.7-1, indicating a 
hard X-ray tail. These features are summarized 
in Table 1. 
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We also show in Fig. 2 that / and y of the stellar- 
mass black hole are the same as those of super- 
massive black hole at low accretion rates, though 
the radiation pressure-dominant part appears at a 
rate ~ 3 times larger than that of AGN. 



Table 1: Two types of solution 



Spectral 
state 


/ 


Range 
of M 


Pressure 
in disk 




Hard 
Soft 


~ 1 
< 0.5 


Always 
> 0.1M E dd 


pg 

^disk 
pr 

^disk 


~ 1 

> 2 



As we pointed out in Sect. 3, there always exists 
gas pressure-dominant solution with Tdisk ~ 10 4 K 
down to the last stable orbit regardless of accre- 
tion rates, as long as a very large fraction of en- 
ergy is channeled to the corona (Similar result 
is obtained in other way by Nakamura & Osaki 
1993 but for an assumed / value), and radiation 
pressure-dominant solutions are found only at high 
accretion rate with a large fraction of accretion en- 
ergy dissipated in the disk. Therefore, there are 
two solutions in the inner region if the mass accre- 
tion rate is high enough: large / corresponds to 
the hard spectrum and small / the soft spectrum. 
This might explain why some soft X-ray transients 
(SXTs), such as GS 2023+338, continuously stay 
in the hard state throughout the outburst, while 
in other SXTs the spectral transition occurs from 
soft to hard during the decay and/or from hard 
to soft during the rise of an outburst (Tanaka & 
Shibazaki 1996). Observations in AGNs also show 
weaker contributions to the hard X-rays for ob- 
jects with higher accretion rate like narrow-line 
Seyfert 1 galaxies which are thought to have large 
M/M (Boiler 2000 and references therein). We 
thus propose that the way of spectral state tran- 
sition is not unique: a) at very high luminosity 
(L ~ LEdd), both the corona and disk are strong, 
which might be time-dependent; b) at moderately 
high luminosity (L > O.lLEdd), the system is usu- 
ally dominated by disk, for which the spectrum is 
soft; but hard state is also possible, in which only 
a little accretion energy is left in the disk (this is 
different from the proposition by Narayan et al. 
(1998) who claim that only the soft state is re- 



alized); c)at low luminosity (L < O.lLEdd), most 
of the accretion energy is transfered to the corona 
and the power-law type hard X-ray spectrum is 
produced. 
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Fig. 1. — Coronal quantities along distance from a 1O 8 M black hole at accretion rate of 0.1MEdd(= ^Edd/c, 
left panel) and 0.5MEdd (right panel). 
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log L/L Edd 

Fig. 2.— The fraction of total luminosity released in the corona (solid curves) and the flux weighted 
Compton y parameter y (dotted curves) in dependence on the mass accretion rate for a 1O 8 M black hole 
and a 1OM black hole. Remarkably, there is no difference between the cases with different black hole masses 
at L < 0.1L Edd . 
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